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Using a three-dimensional (3d) effective field theory and non-perturbative lattice sim- 
ulations, we study the MSSM electroweak phase transition with two dynamical Higgs 
doublets. We first carry out a general analysis of spontaneous CP violation in 3d two 
Higgs doublet models, finding that this part of the parameter space is well separated 
from that corresponding to the physical MSSM. We then choose physical parameter 
values with explicit CP violation and a light right-handed stop, and determine the 
strength of the phase transition. We find a transition somewhat stronger than in 2- 
loop perturbation theory, leading to the conclusion that from the point of view of the 
non-equilibrium constraint, MSSM electroweak baryogenesis can be allowed even for a 
Higgs mass m# ~ 115 GeV. We also find that small values of the mass parameter tua 
(< 120 GeV), which would relax the experimental constraint on ra#, do not weaken 
the transition noticeably for a light enough stop. Finally we determine the properties 
of the phase boundary. 
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1. Introduction 

The observed existence of a baryon asymmetry in our Universe is direct evidence for 
physics beyond the Minimal Standard Model. Indeed, the electroweak theory contains 
anomalous baryon number violation which is rapid enough to be in thermal equilibrium 
at temperatures T > 100 GeV fl], |[ [|, || , so that any pre-existing asymmetry is washed 
out. (Unless there is an asymmetry in baryon minus lepton number, B — L, which 
would also require physics beyond the Minimal Standard Model; for an overview on 
recent scenarios based on this, see [[J). As the Universe then cools down, it turns 
out that there is no electroweak phase transition at all for Higgs masses m# > 72 
GeV [[], |7], , or ra# > 90 GeV in the presence of primordial magnetic fields || : the 
cosmological evolution is just smooth and continuous. Taking the experimental lower 
bound mu ^ HO GeV into account [[l(J (the factual bound is even a few GeV higher 
by now), this would mean that the baryon symmetric high temperature state simply 
freezes in, in contradiction with observation. 

It is quite interesting that already one of the simplest extensions of the Minimal 
Standard Model, the Minimal Supersymmetric Standard Model (MSSM), offers a def- 
inite and experimentally testable way of changing this conclusion. One can uniquely 
identify a bosonic degree of freedom, the right-handed stop, which can be "light" and 
dynamical at the phase transition point without violating experimental constraints 
at zero temperature, and interacts strongly enough with the Higgs to strengthen the 



phase transition significantly [ll]-[]23|. There are also new sources of CP violation 



available which could potentially have a favourable effect pi|-||28||. Many details of the 



electroweak phase transition in this region have recently been studied ||29||-p6|. 

In this paper, we address several issues related to the electroweak phase transition 
in the MSSM. The first is CP violation in the background configuration related to 
the fact that there are two Higgs doublets. The second is the strength of the phase 
transition at an experimentally viable parameter point corresponding to a Higgs mass 
of about 105 GeV, not excluded for the MSSM with smallish tua- (We also consider 
other values close to these, notably van up to ~ 115 GeV.) The third is the structure 
of the phase boundary, or bubble wall, at the physical parameter point. 

In particular, as to the first of these issues, we will pay some attention to the phe- 
nomenon of spontaneous CP violation. Spontaneous CP violation can in principle take 
place in any two Higgs doublet model p7| , but for the physical MSSM parameters it 
cannot be realized at T = |38|, [39j . However, it has been suggested that it might be 



more easily realized at finite temperatures PUL HI] , or even only in the phase boundary 



between the symmetric and broken phases J|2], |43| (in which case it is sometimes called 



"transitional" CP violation). The existence of spontaneous CP violation would mean 
that even small explicit phases can have a large physical effect, and such a situation 
within the phase boundary would conceivably be useful for electroweak baryogene- 



sis 
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The method we use to study all the three questions is the construction of an effective 
3d field theory with the method of dimensional reduction, and its non-perturbative 
analysis with simple lattice simulations. The dimensional reduction step was already 
carried out in |49|], and here we concentrate on the non-perturbative part. 

The plan of this paper is the following. In Sec. |^ we review briefly the form of the 
3d effective field theory for the MSSM, arising after dimensional reduction. In Sec. |3| 
we discuss the lattice implementation of this theory — this section should be skipped 
by those not interested. In Sec. [| we carry out a general analysis of the phase diagram 
of the 3d theory, with particular view on spontaneous CP violation. In Sec. |] we focus 
on a physical choice of parameters and determine the strength of the phase transition, 
while in Sec. |6| we determine the properties of the phase boundary appearing in the 
physical transition, checking also for the possibility of transitional CP violation. We 
summarise and discuss implications in Sec. |7|. 



2. The effective theory 

2.1. The action 

At finite temperatures around the electroweak phase transition, the thermodynamics 
of the MSSM can be represented by a 3d effective field theory containing two SU(2) 
Higgs doublets H 1 ,H 2 and one SU(3) stop triplet U p5| . The action is of the most 
general gauge invariant form, 

£ 3d = \TTGl + {DtU)\DIU)+ml{T)U^U + \ u {U^Uf 

+ ntfUHlHx + l2 U ] UH\H 2 + [ ll2 U^UH\H 2 + H.c. 

+ ^Tr4 + {DfH^^DfH,) + {DfH 2 )\DfH 2 ) 

+ ml{T)H\H x + m\{T)H\E 2 + \m\ 2 {T)E\H 2 + H.c. 

+ \x{E\H x f + X 2 {H ] 2 H 2 f + XsHlHxH^ + X 4 HlH 2 HlH t 



+ 



X 5 (HlH 2 ) 2 + XeHtHiHtHz + X^H^l^ + H.c. 



(2.1^ 



Here Df, Df are the spatial SU(3) and SU(2) covariant derivatives, Gy, F^ the cor- 
responding field strength tensors, and H 2 = ia 2 H 2 . We denote the SU(3) and SU(2) 



gauge couplings by g s ,g w - We have neglected the dynamical effects of the U(l) sub- 
group since they are expected to be small |j0| (although some aspects of the system 
with a dynamical U(l) remain to be understood |^TJ). Note also that compared with 
the MSSM, U denotes the complex conjugate of the original right-handed stop triplet. 
The couplings in Eq. (2.1) can be expressed in terms of the physical parameters of the 



MSSM and the temperature T, as will be specified below. 

For future reference, let us recall that if one is only interested in the strength of the 
phase transition, the effective theory in Eq. ( |2.1|) is even unnecessarily complicated. 
Indeed, it is easy to understand |52|, |TJ|, [15], [16|, [T7|, ^9| (see also [Appendix A| ) that if 
the two Higgs doublet mass matrix is diagonalized, one of the eigendirections is always 
heavy, and can be perturbatively integrated out. This results in an effective theory 
with a single SU(2) Higgs doublet, and the right-handed stop. We studied that effective 



theory with lattice simulations in |19fl . The reason we keep here both Higgs doublets 
is that we measure a number of observables which are specific to the existence of both 
fields, and cannot be addressed with the simpler theory. 

The couplings appearing in the 3d theory in Eq. (|2.1| ) have the dimension GeV, and 
the fields have the dimension GeV 1 ' 2 . For later convenience, we will parameterise the 
couplings by introducing what from the point of view of the 3d theory is an arbitrary 
scale, the temperature T. We then scale all the couplings to a dimensionless form by 
dividing with T, and all the fields into a dimensionless form by dividing with T 1 / 2 : 

2 _ fhd -v _ ^i,3d _ 7i,3d Hj _ -^i,3d U_ _ ^3d /<-, <->x 

9 — rp i 'n — T 1 ' '* — T 1 ' T 1 — Tl/2 ' T — T"l/2 ' \ ) 

Expressed in terms of the newly defined coupling constants and fields, the action goes 
over into 

S = Jd 3 x £ 3 d,ori g -> ^ Jd 3 x £ 3d , (2.3) 

where the dimensionality of the new C^d is GeV 4 as usually in 4d. We shall use this 
formulation throughout the paper, with T taken to be the physical temperature. 

2.2. Approximate physical values of couplings 



Expressions for the parameters in Eq. ( |2.1| ), corresponding to the MSSM, have been 
derived in pT|; for a summary, see Appendix A. 7 there. We cite the precise values used 
in Sees. |5], ^| later on, but let us recall the general orders of magnitude already here. 

It is important to keep in mind a basic difference between the effective theories 
corresponding to the Standard Model and the MSSM. In the former, the physical Higgs 
mass is determined by the effective quartic scalar coupling, A ~ ((7 2 /8)(m 2 7 /m 2 y ), while 
the temperature is determined by the 3d scalar mass parameter, m| d (T) ~ —m 2 H /2 + 



0.3T 2 . In the MSSM, on the contrary, the quartic Higgs couplings are fairly inert, 
~ 5^/8, and are affected by the zero-temperature Higgs mass (i.e., by tan/?) only 
through small radiative corrections. Thus both the physical Higgs spectrum and the 
temperature reside now dominantly in the effective mass parameters. The rough generic 



orders of magnitude for the quartic couplings in the MSSM can therefore be cited J49 
independent of the Higgs mass and temperature: 

Ai ~ 0.07, A 2 ~ 0.13, A 3 ~ 0.08, A 4 0.22, (2.4) 

|A 5 |< 0.002, |A 6 |<0.01, |A 7 |<0.02, (2.5) 

Xu ~ 0.19, g 2 w ~ 0.42, g 2 s ~ 1.1. (2.6) 

In Ai, A3, A4, we have included only the tree-level terms, but in A2 also the dominant 1- 
loop terms proportional to the top Yukawa coupling to the fourth power, hf, which are 
absent in Ai, A 3 , A 4 (we recall that h t « 1.0). In order to get the estimates for A5...A7, 
we have taken into account that when the right-handed stop is light, the squark mixing 
parameters cannot be too large compared with the left-handed stop mass, because of 
the stability of the theory (see below). 

As to the three couplings 7$ in Eq. ( |2.1| ), they can be reparameterised in terms of 
the three couplings h t , A t , \x as 

71 = -h\\ii\\ l2 = h\(\-\A t \ 2 ), 7i2 = /* 2 4*£*, (2.7) 

where A t , jjl are the mixing parameters in the 3rd generation squark mass matrix, rriQ 
is the corresponding left-handed squark mass parameter, and A t « A t /mQ, fi « n/rriQ. 
There are of course radiative corrections to these relations, but we can also view them 
as a more general reparameterization, since At, jl are essentially free parameters. We 
mostly assume A t , fi< 0.3, again in order not to destabilize the theory (see, e.g., fl49"|). 
and also since large values tend to weaken the phase transition (see, e.g., JTJ], ^, |20|). 
Consider finally the mass parameters. Working in the limit rhfj = —mfj -C (27rT) 2 ^C 
ttiq, where mfj is the right-handed stop mass parameter, we have at leading order 

ml{T) « -m^+^ + ^-^dAp+l/il 2 ))^, (2.8) 



m\{T) « ^ml-^(m\ + ml)cos2(3 + (^gl-hM 2 y 2 , (2.9) 

m\{T) « imU^^ + ^^^ + g^ + ^-i^li.l^T 2 , (2.10) 

m\ 2 (T) « -\m\mi2l3 +h 2 t A* t fi*T 2 . (2.11) 



Here m^,tan/3 are the usual MSSM input parameters. We cite these expressions 
because they lead to some generic properties relevant for our discussion below, in 
particular that the trace of the two Higgs doublet mass matrix, m\(T) + m\(T), is 
always positive in the region relevant for us. 



3. Lattice formulation 

For future reference, we recall next briefly how the theory in Eq. ( |2.1| ) can be discretized. 

3.1. Lattice action 

We discretize the action in Eq. fl2.1|) in the standard way. The finite lattice spacing 
a enters as f3 w = 4/(ag^T), (3 S = 6/(ag 2 T), and the lattice volume is denoted by 
N 1 N 2 N 3 . The gauge field terms (l/2)TrG?. + (l/2)TrF? are treated with the usual 
Wilson formulation, as in [H5|. The scalar fields are rescaled into a dimensionless form 



by Hi^Hi = Hi/T, U ^U = U/T, and then, e.g., 

S = ^Jd'xliDtUViD^+mUT^U + XuiWUf 
]T[-2aT]TRe^A + i 

x i 

+ (6aT + (aT) 3 (m 2 UtB /T 2 ))u$ x + A a (aT) 3 (^C/ x ) 2 ]. (3.1) 

Here U^ i is the SU(3) link matrix at point x in direction i, and the bare lattice mass 
parameter mfj B is given in [Appendix B . 



Let us note that by computing mfj B we have fully renormalized the theory [|53" |, but 



we have not implemented 0(a) improvement |54| here. This is because according to 



Sec. B.4 of the latter ref. in [Q, one would need /3 W > 30 to be comfortably in the 
regime where improvement works, and we are not able to go that close to the continuum 
limit, with the lattice sizes we can manage in practice. Thus we expect a more general 
lattice spacing dependence (and use correspondingly a more general fitting ansatz for 
the continuum extrapolation). 

3.2. Update algorithm 

The lattice simulations of the theory in Eq. (|2.1|) are quite demanding, due to the 
large range of mass scales near the first order transition temperature (we recall that 
many mass scales have already been removed by the analytic dimensional reduction 
computation, but a number of dynamical scales are still left over, particularly since we 



want to keep both Higgs doublets in the effective theory). The Monte Carlo update 
algorithms employed are nevertheless qualitatively similar to the ones used in ref. |19| 
for simulating MSSM with only one SU(2) Higgs doublet. 

Our update algorithm consists of a mixture of overrelaxation, heat bath and Metropo- 
lis updates for both the Higgs and gauge fields. For all of the three Higgses (Hi, H 2 , U), 



we use the "Cartesian overrelaxation" presented in |55], [n|. The overrelaxation update 
is much more efficient in evolving the fields than the diffusive heat bath and Metropolis 
updates; however, in order to ensure ergodicity one has to mix diffusive update steps 
with overrelaxation. We use a compound update cycle which consists of 5 overrelax- 
ation sweeps through the lattice, followed by one heat bath/Metropolis sweep. For 
details, we refer to Sec. 6 of |B?|. 



For the simulations to be reliable and to allow for an extrapolation both to the 
infinite volume and continuum limits, the lattice spacing a and the lattice size L = Na 
should satisfy two inequalities, 

0- < £min; £max < Na , (3.2) 

where £ min and £ max are the smallest and the largest physical correlation lengths of the 
system. In the present system the mass scales near the transition can vary by roughly 



an order of magnitude (see Sec. |5.4| ). This makes it necessary to use relatively large 
lattice sizes N, but even then it is not easy to satisfy the inequalities in Eq. ( |3.2| ) with 
wide margins. This emphasizes the importance of a careful check of both the infinite 
volume and the continuum extrapolations. 

The fact that the transition is strongly of the first order, makes the inequalities in 
Eq. ( |3.2|) somewhat easier to satisfy than for instance in the case of a second order 
transition, since £ max remains finite. However, the transition is now so strong that the 
system does not spontaneously tunnel from one metastable phase to another, especially 
in large volumes. On the other hand, probing the whole tunnelling phase space region is 
required in order to determine precisely the critical temperature, the order parameter 
discontinuities, and the surface tension. Thus, to allow for frequent tunnellings, we 
use the multicanonical simulation method with automatic weight function calculation, 
discussed in |T^ . 



3.3. Observables 

The simplest observables to be considered are the gauge invariant bilinears 

1 = H\H x , 2 = H ] 2 H 2 , R = *RbH\H 2 , I = lmH\H 2 . (3.3) 




The operator I is odd under complex conjugation, and constitutes thus an order pa- 
rameter for spontaneous C violation (see below). In addition to the global expectation 
values of these operators, we also consider the corresponding correlation lengths, ob- 
tained from the 2-point functions. 

In perturbative analyses, one often uses the SU(2) and U(l) symmetries of the theory 
to parameterise the two Higgs doublets for instance as 

H 2 = ^( COSeel '), (3.4) 

y/2 \ sin9 J ' V ^ 

and tan/3 = v-ijv\. When used beyond tree-level and in connection with, say, some 
covariant gauge, the values of the quantities v%, v 2 , 9, <fi in the broken phase are gauge 
dependent, and thus there is no unique relation to the values of the operators in 
Eq. (|3.3| ). However, at the phase transition point (T = T c ) we may define a gauge and 
scale independent generalization of the perturbative parameters for instance by 

\v\ ee (HjH t ) - (HjH t ) . , tan 2 (5 = %, (3.5) 

2 symmetric ij-? 

where i = 1, 2 and the latter expectation value is taken in the homogeneous "sym- 
metric" high-temperature phase. We could also write = arctan( I/R), or <p = 
arctan[(7 — /symmetric )/(R — -^symmetric)], but such ratios are numerically very unsta- 
ble, and we do not use them here. Rather it is the perturbative values of Vi,V2,9,(f> 
which should be converted into gauge invariant observables such as those in Eq. (|3T 

3.4. Mean field estimates 

Finally, let us note that due to the relatively heavy cost of simulating the action in 
Eq. (|2.1|) on the lattice, many of the preliminary parameter scans have been made using 
relatively small lattice sizes. To at least partly account for the finite volume effects 
in the comparison with perturbation theory, we transform the perturbative results for 
the quantities in Eq. ( |3.4j ) into finite volume "mean field" estimates for the operators 
in Eq. ( |3TB| ) as follows. 

A mean field estimate can be obtained by taking the fluctuations into account in 
the effective potential, and performing then the integral over the zero-modes of the 
fields, parameterised by v\,V2,9,cj) as in Eq. (|3.4|) . In addition, we take into account 
the fluctuations of U, and parameterise U = (l/\/2)(x, 0, 0) T . The integration measure 
goes over into 

dRxdE^dU = C dv dx df3 d9d^v 7 x 5 cos 3 (3 sin 3 Pcos9sin9, (3.6) 



where C is a constant. The action can be written as 

S«^Kff(^AM,x), (3-7) 

where V = N\ N2N3C1 3 is the volume, and the mean field estimates are then obtained 
from 

(O t ) « Z- 1 fdH 1 dH 2 dUO l ex V (-S) 1 (3.8) 

where the 0{ are written using Eq. 



4. Spontaneous CP violation, and transitional 

We now move to the first physics topic. In this section we consider the case of no 
explicit CP violation (i.e., all the parameters in Eq. ( |2.1| ) are assumed real), in order 
to carry out into completion the analysis outlined in p9| . The motivation is that if 
spontaneous CP violation would exist in the system where CP phases are put to zero, 
then one could expect large physical effects once small explicit phases are turned on. 

4.1. Symmetries and phases 

Let us start by reviewing briefly the overall setup. In the space of general couplings, 



the theory in Eq. ( [2.1|) has several "phases". First of all, there are the usual phases 
related to "gauge symmetry breaking": the SU(3) gauge symmetry can be broken by 
an expectation value of U, and the SU(2) symmetry by expectation values of Hi,H2- 
In the physical case the parameters had better be such that one does not end up in 
the phase with a broken SU(3) symmetry, since one would be stuck there forever [^2] . 



In addition to the gauge symmetries, there are also global symmetries in the system 
of Eq. ( |2.1| ). There is one continuous U(l) symmetry corresponding to the hypercharge, 
and then there are the usual discrete symmetries C, P. The time translation symmetry 
T is not directly visible any more in the 3d effective theory, but it dictates what kind 
of operators can arise in the dimensional reduction step |56] . 



As to the continuous U(l) symmetry, let us first recall how things are in the single 
Higgs doublet SU(2) theory, C = (DfH)\DfH) + ... . In this case the global U(l) 
symmetry is H — > exp(ia)H. This global symmetry cannot get broken, however, since 
any configurations H, exp(ia)H can be SU(2) gauge transformed to the same config- 
uration (e.g. to the unitary gauge). Thus the system is U(l) invariant independent of 
the expectation value of H, and there is always a massless excitation in the gauged 
version of the theory. 
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Things are different if there are two Higgs doublets, C = Y J j=i,2{DfHj) J[ {DfHj) + ... . 
Taking into account the general form of the two Higgs doublet potential, Eq. (|2.1|), there 
again remains a global symmetry H\ — > exp(ia)Hi, H 2 — > exp(—ia)H 2 . However, now 
this symmetry can get broken: if Hi,H 2 are not proportional to each other, it is 
not possible to unwind simultaneously the angle from exp(ia)Hi,exp(—ia)H 2 by an 
SU(2) gauge transformation. For such expectation values of H\, H 2 , physics is not U(l) 
invariant. In the gauged version of the theory, the photon becomes massive. In terms 
of the parametrization in Eq. ( |3.4| ), the breaking of the U(l) symmetry corresponds to 
I sin 6*| > 0, or I cos#| < 1. 

As to the discrete symmetries, for real parameters this theory is even under both of 
the discrete symmetries C, P. The C transformation corresponds to 

Hi - H*. (4.1) 

While parity is not expected to be spontaneously broken in this theory []57|, Q, the 



C symmetry can be, thus violating also CP ||37|| . The breaking of C is signalled by 
a non-vanishing expectation value of the local gauge invariant order parameter I in 
Eq. ( |3.3j ). In terms of the parametrization in Eq. ( |3.4j ), C symmetry corresponds to 
the invariance of the theory under — > —(f), and the breaking of C is signalled by 
\cos(p I < 1. 

4.2. The phase diagram in 1-loop perturbation theory 

We now want to find out the parameter values for which the (global) phases discussed 
in the previous section are realized. We first do this in perturbation theory, using 
the (gauge specific) variables in Eq. Q3.4p . We work in the Landau gauge. We al- 
ways assume rriu(T) > 0, in order to avoid the dangerous charge and colour breaking 
minimum 



22|. For real A t , fi, the 1-loop effective potential is then 



V{vi,V2,6,4>) = -m\(T)v\ + -m\{T)v \ + ml 2 (T)viv 2 cos 8 cos 

+ t-Vi + -A 2 ^2 + jhvlvi + j^v \v\ cos 2 9 

1 1 

+ -A 5 v 2 V2 cos 2 9 cos 2<p + -{\v\ + A 7 f 2 )fif 2 cos 9 cos <p 

~ TT-aKvl + vl)! -—(A + Bcos9cos^ --— E(m^)l. (4.2) 

107T Z7T 1/7T ~[ 

Here m 2 Si are the real eigenvalues of the 8x8 scalar mass matrix, obtained after making 
a shift according to Eq. (|3.4|) in Eq. (|2.1|), and 



A = ml(T) + hi\vl - \h\(\li\ 2 v\ + \A t \ 2 v 2 2 ), B = h 2 A t ^ Vl v 2 . (4.3) 
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We now note that the dominant 1-loop effects in this effective potential are the terms 
from the vector bosons and from the stops, while the 1-loop scalar effects from m 2 Si on 
the last line in Eq. (|4.2p are small. This is because the scalar self-couplings A» are never 
large in the MSSM, being <g^/8 (see above). Furthermore, we may set A t = jl = 
in Eq. ( |4.3|) for the moment, which allows for a simple analytic treatment. We then 
work out a complete parametrization for the part of the parameter space leading to 
spontaneous C violation. The same could be done for the phase with broken U(l), 
but as it turns out to lie even farther away from the MSSM, we do not elaborate on 
it here. The details of the discussion concerning the C violating phase are presented 
in [Appendix C| , and we address now the results only. 



To be explicit, we fix the couplings A1...A4 to the values given in Eq. (|2.4|). We then 
make a full scan of the remaining parameter space according to Eqs. ( |U.10D -( |(J.12| ), 
(|C.24|) -( p.26|) , without any additional restrictions on A5, Ae, A7, rn\{T) 1 rn\{T),rn\ 2 {T). 



As to the expectation value v/T, we take the realistic MSSM into account by recalling 
that there one does not get values as large as v/T ~ 3 at the electroweak phase 
transition, and in any case for such vevs dimensional reduction and the construction of 
the effective theory in Eq. ( |2.1| ) start to lose their accuracy. Thus we assume v/T<3. 
The results of the scan are shown in Fig. [I]. We show the projections of the parameter 
space onto different axes, in most cases with \m\(T) + m\(T)\/T 2 as the x-axis. We 
have only shown the part [m\(T) + m 2 (T)]/T 2 > 0, as this is the case relevant for the 
MSSM, see Eqs. ([TSp, ( g7TU|) and the discussion below them. 



Let us make a few observations on the results. First of all, for A1...A4 as given, the 
"bounded from below" constraint in Eq. (|C.9|) forbids values A 5 > 0.026, and the fact 
that we only allow values v/T < 3 leads to a further restriction visible in Fig. [j], cf. 
Eq. (|C.26|) . (Furthermore, note from Eq. fl2.5|) that only small values are produced by 
dimensional reduction.) Second, Mf 2 = rn\ 2 + (1/2) (\v 2 + A 7 t>|) is very small; this is 
in order to satisfy Eq. ( |U.2|) , for the small values of A5 that arise. Finally, the x-axis of 
the figures, [m\{T) + m\{T)]/T 2 1 is also restricted to quite small values. This is again 
ultimately due to the smallness of A5, as it was argued in p9[ that one has to satisfy 



m\{T)+m 2 2 {T)<2\,v 2 . (4.4) 

This is strictly speaking quantitatively true only at tree-level, but we can observe from 
Fig. H that there is no order of magnitude change due to 1-loop effects. 

As we discuss in Appendix |C.6| , the inclusions of 1-loop effects from the scalar self- 



couplings (contained in m 2 Si ) and from allowing A t , fi 7^ in Eq. Q4.3|) do not change 
these results in an essential way. Let us reiterate that the Aj's are always small in 
the MSSM, so that one practically never enters the region relevant for the Standard 
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Figure 1: Projections of the region of spontaneous C violation (| cos 0| < 1) onto 
various planes. Here M\ 2 = m\ 2 {T) + (l/2)(\§v\ + A7t>§), M = \/2mjj(T)/h t , cosa 6 = 
A 6 /(2v^A;), cosa 7 = A 7 /(2v^A7), cf. Eqs. Q, (|Cl3|) , (|02lD , (|025D . Thin 
contours (equally spaced, arbitrary normalization) indicate roughly the thickness of 
the allowed region of the parameter space in the orthogonal directions. 
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Model and many other systems, where scalar fluctuations related to Aj's change the 
predictions of perturbation theory in a qualitative way (see, e.g., || |58|, |59|). 



4.3. The non-perturbative phase diagram 

We now wish to explicitly check how the perturbative predictions above are changed 
by non-perturbative effects. The general experience from 3d gauge+Higgs systems is 
that the properties of phase transitions are badly described by perturbation theory 
if the transition is weak so that the smallest mass scale m m j n appearing is small, 
g^T / (7rm m i n ) > 1 J55J . This happens typically for large scalar self-couplings. On 
the other hand, the changes in the locations of the phase transitions are expected 
to be small even in that case: in the critical temperature T c , non-perturbative ef- 
fects arise only at next-to-next-to-leading order |Kj, so that parametrically 5[m\{T c ) + 
m 2 (T c )}/T 2 ~ ^g^,/(47r) 2 . We would now like to verify how well this is true numerically 
and, in particular, whether C violation could be more favoured than was perturbatively 
estimated. (For example, could it take place at somewhat larger \m\(T) + m\{T)]/T 2 l) 
In order to carry out this check, we shall increase [m\(T) + m\(T)]/T 2 starting 
from the phase where C is broken, crossing the phase transition. We parameterise the 
starting point as in Eqs. (|C.10|) -( |C.12|) , and add then a new dimensionless parameter 

ato m ?(r)/r 2 ,m|(r)/r 2 : 

nm^ap +v , i=1 , 2 . (4 . 5) 

For simplicity we shall also denote what used to be v/T in Eqs. ( |C.10| )- (fC.12| ) by 
v. Increasing y corresponds to increasing the temperature in the 4d language, while 
v, tan/5, cos0 are used to parameterise rn\{T)/T 2 ,in\{T)/T 2 1 rn\ 2 {T)/T' 2 ' at the refer- 
ence point y = 0. 

The parameter space of the theory is quite large, so we have to make some choices. 
In the following, we choose a relatively small value A5 = 0.001 because this is natural 
from the point of view of dimensional reduction in the MSSM (cf. Eq. Q2.4J) ), and since 
a small value leads more easily to relatively small v/T, not much larger than unity (cf. 
Eq. (|U.26 )), which is also what we expect around the electroweak phase transition. We 
do not expect our results to change qualitatively with variations of A5. 

We will also fix cos0 = 0.5, which does not have a significant effect. Moreover, we 
choose here a 6 = tt/2, 07 = n/2, setting A 6 , A 7 to zero (according to Eqs. ( |C24| ), ( |U.25|) , 



A 6 = 2VA1A5 cos ae, A7 = 2v / A2A5Cosa7). If one would take other values, one needs to 
have an anti-correlation in the signs of A6, A7, to get Mf 2 = m i2(^) + \^g v i + \^"t v \ 
close to zero; see Fig. [I]. However, we again do not expect new qualitative effects 
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\ 5 =0.001, X=X 7 =0.0, cos^)=0.5, M/T=0.5 
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Figure 2: The 1-loop phase diagram for the transition from the C violating phase to 
the usual broken and symmetric phases as y is varied, for fixed A5 = 0.001, A6 = A7 = 
(« 6 = a 7 = tt/2), cos0 = 0.5, M/T = 0.5 ( mu {T)/T = 0.35), A t = fi = 0. The 
transition from the C violating phase can lead either directly to the symmetric phase 



'"s." x 



or to the usual broken phase ("b."), and it can be either of 1st or 2nd order. 



from relaxing this assumption. As to the stop sector, we somewhat arbitrarily choose 
ml(T)/T 2 « 0.13, but vary A t , jl within \A t \, |£| <0.3. We denote M = \f2m u \T)/h t 
(cf. Eq. ( |Oll )) so that M/T = 0.5. 

These choices fix the quartic couplings of the theory as well as the stop sector, 
but still leave the diagonal entries of the SU(2) scalar mass matrix, parameterised by 
v, tan/? in Eqs. ( |C.10| ), ( |C.11| ), open. The perturbative phase diagram in this space, 
based on the full 1-loop effective potential in Eq. (|4.2|), is shown in Fig. 0. We have 



chosen a number of points from this diagram for further non-perturbative study, such 
that all different qualitative types of transitions are represented. 

In Fig. [| we show the mean field estimates for the behaviour of the operators in 
Eq. ( |3.3|) for a few representative choices of parameters. These estimates follow from 
Eq. (|3.8|) , with the 1-loop effective potential from Eq. ( |4.2j ), supplemented by a non- 
vanishing value of x as needed in Eq. ( |3.8|) , whereby the contribution on the last line in 
Eq. (|4.2j ) goes over into a sum over the eigenvalues of a 9x9 matrix, with the radial U 
direction coupled to the SU(2) scalars. (The Goldstone modes of U appear separately 



as, e.g., m 



00 



The mean field estimates can be compared with the lattice results, shown also in 
Fig. |3| for the same parameter values. We observe that, as expected, the behaviours 
are quite close to each other. We then estimate that compared with Fig. [l], the largest 
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Figure 3: Top: Mean field estimates for the ^-dependence of different expectation 
values, corresponding to the physical volume obtained with j3 w — 8, NiN 2 N 3 = 20 3 . 
For / the absolute value of the volume average is taken. Bottom: Lattice results at 
P w = 8. For the first two sets the transition is strongly of the first order, and a small 
volume has been used; in spite of this, some y- values allow for separate measurements 
in two different metastable phases. For the middle set, the transition goes on the lattice 
directly to the C conserving symmetric phase, rather than to the C conserving broken 
phase seen in the perturbative plot. 



value allowing for spontaneous C violation can change at most by 
where the upper bound is quite conservative. 



(4.6) 
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4.4. Implications for the physical MSSM 

Let us now compare Fig. [l] supplemented by Eq. Q4.6Q , with the part of the parameter 
space allowed by the MSSM. One could make many comparisons, but we focus here on 
m\ 2 (T)/T 2 and, in particular, [mj{T) + m 2 2 {T)]/T 2 , cf. Eq. fl4j ). 

First of all, we recall from Fig. [1] that one needs small values of M\ 2 jT 2 and therefore, 
due to the smallness of A 6 ,A 7 , small values of m\ 2 {T)/T 2 (cf. Eq. (|C.1| )). It can be 
observed from Eq. ( j2.11| ) that in the MSSM this is easier to satisfy at finite temperatures 



than at zero temperature, due to a temperature correction which can cancel the T = 
part in m\ 2 (T) for A t fi > 0. Furthermore, even if the experimental lower limit on m^ 
appears to be rather high, ~ 90 GeV |10|| , one can at least partially compensate for 



this by taking a large tan/3, since only the combination m^tan/3/(l + tan 2 /?) appears 
in Eq. (|Tl|). 

However, the constraint on [mf(T) + m 2 (T)]/T 2 works in the opposite direction. 
To get spontaneous C violation m\{T) + m 2 {T) should be small (Fig. |l|), its order 
of magnitude given by A5V 2 (Eq. ( |4.4| )). Finite temperature does not help with this 
constraint at all: v 2 gets smaller, while m\(T) + m 2 (T) gets larger. 

To be more precise, we obtain using Eqs. (|IO|), ( |2. 1Q ) that in the MSSM, 

m\(T) + m 2 2 (T) &m 2 A + 0.5T 2 + 0.25T 2 (1 - |i t | 2 - \fi\ 2 ). (4.7) 

Let us now reiterate that in the limit of large mq and small mfj that we are working 
at, 1 — \A t \ 2 — \fi\ 2 should in general be positive for the theory to be consistent from 
the point of view of boundedness and electroweak vacuum stability |49| . Taking into 



account also the experimental lower limit on the CP odd Higgs mass parameter vha ^ 90 



GeV 0, we then get that in the MSSM, 

K(T)+m 2 (T)]/T 2 > 1.3. (4.8) 

This holds for all temperatures below 100 GeV, i.e., also in the broken electroweak 
Higgs phase. 

This result can be contrasted with Fig. |l|. We observe that there is no overlap, since 
[mf(T) +m 2 (T)]/T 2 < 0.5 always. A non-perturbative change of the order in Eq. ( |4.6|) 
clearly cannot bridge the gap. 

It can be observed from Fig. |l]that increasing mu{T) seems to allow for larger values 
of [m\{T) + m 2 (T)]/T 2 . However, this effect is not enough to change our conclusions. 
In fact, for large mjj(T) the field U can be integrated out, as we review in [Appendix D , 



and the result is a theory of the form in Eq. ( |2.1| ) but without U . In this theory, there 
is an upper limit on [m\{T) + m 2 {T)]/T 2 leading to spontaneous C violation as well, 
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Figure 4: The tadpole graph leading to Eq. (|4.10|) . The effective local quartic vertex 
shown, the second line in Eq. ( |2.1| ) with 7j's from Eq. Q2.7|), is a good approximation 
as long as mfj -C thq. 



numerically [m\{T) + m 2 2 (T)}/T 2 < 0.4 for v/T< 3.0 0. The effect observed in Fig. |T| 
is equivalent to the fact seen in Eqs. ( P.l|) , (|D.2| ) that the effective [m\{T) + m 2 ] (T)}/T 2 
after integrating out U tends to decrease by increasing mu(T). More precisely, the last 
term in Eq. (|47|), (/i 2 T 2 /4)(l - \A t \ 2 - |/a| 2 ), is multiplied by the factor 

/ = 1 - ^Si. (4.9) 

TV 1 

However, this decrease does not continue forever, but the formula breaks down for large 
enough mjj(T) when the high temperature expansion is no longer applicable, and / 
then goes over into the tadpole integral shown in Fig. [| (provided that still mfj -^ ttlq), 



, 6 f°° , x 2 
J -»• -» / ac 



7T 



Va; 2 + ?/ 2 exp(Va; 2 + y 2 ) — 1 



(4.10) 

y=m u /T 



This is always positive, so that [m\{T) + m\{T)]/T 2 does not decrease below m 2 A /T 2 + 
0.5 > 1.3 even if mjj{T) > T, and our previous conclusion continues to hold. 

In summary, spontaneous CP violation seems to be excluded in the MSSM both at 
T = and at finite temperatures around the electroweak phase transition. 

4.5. Transitional CP violation 

Let us finally come to the issue of "transitional CP violation". There have been sug- 
gestions that even if not in the symmetric or broken phase, spontaneous C viola- 
tion could take place within the phase boundary between the symmetric and broken 
phases [JO, ^3[. However, these suggestions could not be confirmed by later analyses 



for physical parameter values (particularly realistic tha) [j3U| j. Basically, the point is 
that as our discussion above indicated, spontaneous C violation is always more likely 
at large vevs, cf. Eqs. ( |C.2|) , (|4.4| ). Thus it seems unlikely that C would be violated in 
the phase boundary, if it is not violated in the broken phase. Below, we will inspect 
the issue numerically at a physical phase boundary corresponding to the MSSM, and 
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find the same conclusion. Let us note that, on the contrary, transitional C violation 



could take place in, say, the NMSSM pO, |29 



5. The strength of the physical phase transition 

5.1. Parameter values 

We now move to the electroweak phase transition in the physical MSSM. Let us start by 
choosing parameter values. We have previously carried out simulations corresponding 
to tan/3 = 3,rriH = 95 GeV, a large m^ > 200 GeV, and a light right-handed stop J05). 
We found a transition which was somewhat stronger than in 2-loop perturbation theory 
and certainly strong enough for baryogenesis. Recently, 4d finite temperature lattice 
simulations have also been carried out in a scalar theory with MSSM type couplings, 
at ran ~ 45 GeV |62fl . There the transition is very strong, and it was found to agree 
well with perturbation theory. 

We now want to take a larger tan/3 = 12 than in |Tj|, corresponding to mn ~ 105 
GeV for a left-handed squark mass parameter niq « 1 TeV and a light right-handed 
stop. We also take a smaller vha ~ 120 GeV for two reasons. First, because then the 
experimental Higgs mass lower bound is somewhat relaxed [pTD| . (Recently it has been 
suggested that the experimental Higgs mass lower bound may be further relaxed for 



such vriA because of large explicit CP violation in the system, as we will have |p3|-||66||). 
Second, because a smaller txia makes the heavy Higgs doublet eigendirection somewhat 
more dynamical, since effects related to it are suppressed by ~ g 2 w T /{m 2 A + 0.5T 2 )2, 
see |Appendix A| . This means that the situation could be favourable for "dynamical" 
CP violation ||9| (i.e., a somewhat non-trivial profile for CP odd observables within 
the phase boundary, even if not actual spontaneous CP violation), as well as for having 
a non-trivial tan /3-profile |67], ^0[ |3(J, which might affect the actual baryon asymmetry 



produced 0-^| 

We also introduce a non- vanishing squark mixing parameter A t w 200 GeV, as well as 
gaugino and Higgsino mass parameters M2, fi, with M2 ~ |/i| ~ 200 GeV. The strength 
of the phase transition depends little on these parameters [pl]1-[p3|. In addition, to 



observe the CP violating effects more clearly, we choose a maximal explicit phase for 
the /i-parameter, \x = i\fi\. The first and second generation squarks and sleptons are 
assumed heavy, whereby even such a large phase is not in conflict with electric dipole 
moment constraints []68|-|[72|. 

Finally and most importantly, we take a negative right-handed squark mass param- 
eter, rhfj = — mfj > 0. Most of the time we work at rh\j = 65 GeV (see below). To 
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summarise, we thus have 

(£«0.42, gjal.1, ^ « 1.0, (5.1) 

m Q « 1 TeV, my « 65 GeV, m A « 120 GeV, tan/3 « 12, (5.2) 

At » |/i| « M 2 » 200 GeV, ju««|/*|, (5.3) 

where on the first row the couplings are assumed to be evaluated at a scale ~ 2-kT . 
These parameters correspond to a lightest physical Higgs mass of about 105 GeV, and 
a lightest stop mass of about 155 GeV, with an uncertainty of a few GeV. 

Applying then the formulas in Appendix A. 7 of ||9| (and fixing T ~ 95 GeV, as 
suggested by 1-loop perturbation theory, inside logarithms and elsewhere where its 
effects are subdominant, in order to simplify the results), we obtain the following 
effective couplings for the theory in Eq. Q2.1|) : 



m\{T) « 18380 GeV 2 + 0.1218T 2 , (5.4) 

m\(T) « -3980 GeV 2 + 0.6218T 2 , (5.5) 

m\ 2 (T) « (-1190 - zlOO) GeV 2 + i0.0030T 2 , (5.6) 

mKT) « -4225 GeV 2 + 0.8534T 2 , (5.7) 

7i « -0.044, 72 « 0.96, 712 « -i0.037, A^ « 0.197, (5.8) 

Ax « 0.0652, A 2 « 0.1188, A 3 « 0.0673, A 4 « -0.1948, (5.9) 

A 5 ^ -0.00019, A 6 ^i0.0017, A 7 « i0.0022. (5.10) 

We should stress that these numbers have of course some perturbative errors, but this is 
not essential for our main statements. Indeed, we will compare 3d perturbation theory 
with 3d lattice simulations, and precisely the same parameters in Eqs. ( |5.4|) -( [5~T0| ) are 
chosen in both cases. This will allow us to unambiguously find out whether there are 
non-perturbative effects in the system. Such non-perturbative effects will then remain 
very similar even if the 4d parameter values are changed slightly, or if the reduction 
computation leading to Eqs. (|5T4])-( |5TT0"| ) is carried out more precisely. 



Finally, let us mention a technical point. We treat the mass parameters in Eqs. 
( |5.7| ) as those at the scale p, = T inside the 3d theory (to be more precise: we choose 
the A-parameters discussed in Appendix |A.4j to be A ~ LOT). In order to remove 
the ambiguity from this choice, a complete 2-loop dimensional reduction computation 
would be needed for the mass parameters [J52]] . Unfortunately, such computations have 



been carried out only in special cases. One is the Standard Model [|5^| , where it turns 



out that numerically A ~ 7T. In |1| it was argued that this should be expected also for 



the MSSM. An explicit computation was then carried out in [21]] for a small tuq ~ 300 
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Figure 5: The perturbative 2-loop phase diagram for niq = 1 TeV. The values of v/T 
at the transition point are also shown. Lattice results are displayed at rhu = 65 GeV 
(square), and at the triple point (triangle). The lattice triple point errorbars include 
only statistical errors (see the text), and are thus an underestimate. 



GeV, and the actual scale was found to be of order ~ IT for the diagonalized Higgs 
mass parameter (see [Appendix A| for the definition), ~ IT for the stop mass parameter. 
However, the scales depend on the other parameters of the theory. While the way to 
carry out the computation for large ttiq, the case relevant here, has also recently been 
clarified f73fl , explicit results for the full MSSM are still missing, so we cannot simply 
take over old values. 

Fortunately, it turns out that this ambiguity is of very little significance for our 
results. Indeed, we have tested the effect of changing from A ~ T to A ~ IT with 
2-loop perturbation theory inside the 3d theory. We find that the critical temperature, 
as well as the value of rhu corresponding to the "triple point" (see Fig. ||) change 
by a few GeV, but apart from this shift the values of v/T at the transition remain 
almost the same. Thus the ambiguity is completely inconsequential, if we normalise 
our parameter values with respect to the triple point. 



5.2. Phase transition in perturbation theory 



For the parameter values in Eqs. ( p.4| )-( |5TT0| ), the theory in Eq. ( |2.1| ) has a first order 
electroweak phase transition. Let us first discuss its properties in perturbation theory. 
In the following, we use the Landau gauge and the scale parameter p, — T, as is 
conventional in the literature. 
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Figure 6: The properties of the phase transition in 1-loop and 2-loop perturbation 
theory, as well as on the lattice. Left: the perturbative effective potentials at the cor- 
responding critical temperatures. Right: the vacuum expectation values as a function 
of T. The lattice value refers to w/T| latticc = (2AJ2 i HJH i )i, cf. Eqs. flOD, 



As we mentioned in Sec. [2.1| , for studying the strength of the phase transition the 
theory in Eq. ( |2.1| ) can be simplified by integrating out a linear combination of the two 
Higgs doublets. This is not only a convenience but also a way of increasing the accuracy 
of perturbation theory: large effects related to a heavy excitation get resummed. We 
discuss the details of the integration out in [Appendix A[ After the integration out, 
we can directly employ the 2-loop effective potential computed in [|19| . We may note 
that at 1-loop level we have also the effective potential in the full theory available, see 
Eq. (|4.2|), and in practice we find quite similar results as by using the diagonalized 
effective theory (to 1-loop accuracy). 

In Fig. |5] we show the phase diagram as a function of mjj,T. The three familiar 
types of transitions [T7| are well visible. However, with our parameters, we observe 
that in fact a 2-stage transition as proposed in [jT7[ is excluded, unlike with the pa- 
rameter choice which we employed in |lj|. Instead one has to worry about whether 



the metastability of the broken Higgs phase is sufficiently strong [ 731 , if rhu is very 



close to the triple point. However, we can here apply the logic of p2| in the opposite 
direction, and state that one would never tunnel into the broken U direction even if 
possible in principle. (In the case of [22], the dashed line in Fig. [5] was tilted in the 
other direction, and the statement was that if one ends up in the phase with broken U, 
one can never get from there to the usual electroweak phase, even if that would be the 
thermodynamically stable phase.) Still, values rhu % 70 GeV would seem safest. We 
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choose here fhu = 65 GeV for further study. 

As discussed in Sec. |3]3|, the strength of the phase transition is in perturbation 
theory usually addressed in terms of v/T rather than H\Hi + H 2 H 2 , or Wh in the 
diagonalized theory of Appendix A . The relation is then obtained from Eq. ( |3.4|) or, if 
one converts from lattice to perturbation theory, from Eq. ( 13.51) . We show the 1-loop 



and 2-loop results for v/T in Fig. |6|. We observe the familiar feature that the transition 



is significantly stronger at 2-loop level |13|-pi|, which is related to the fact that the 
critical temperature T c is lower. 

Finally, let us recall that because a certain combination of the Higgs doublets is 
heavy, the transition takes essentially place in one "light" direction in the field space 
much like in the Standard Model (see also Fig. [T3| below), and the sphaleron energy is 
really determined by the value oi HlH 1 +H 2 H 2 in the broken phase. Note, in particular, 
that even though possible in principle |75| , we would not expect the sphaleron bounds 
to be modified by the CP violation apparent in the couplings in Eqs. ( E5.4|) -( l5Jl)| ), 
because the Higgs direction orthogonal to the light eigenmode h is indeed so heavy 



that effects related to it are strongly suppressed; see [Appendix A . 



5.3. Lattice simulations 



Pv, 


volumes 


8 


12 3 


16 3 


24 3 


12 


16 3 


24 3 


32 3 


16 


24 3 


32 3 




24 


32 3 


40 3 


56 3 



Table 1: Lattice spacings (3 W = 4/(p^aT) and lattice volumes, used for simulations at 
the transition temperature. All the simulations here are multicanonical. 



In order to inspect the reliability of the perturbative estimates discussed above, we 
have carried out lattice simulations at fhu = 65 GeV. First, a series of simulations was 
performed in order to determine the transition temperature T c , the value of v/T in 
the broken phase, the latent heat and the surface tension. The lattice sizes and lattice 
spacings are shown in Table [I]. For each lattice listed we performed 50 000 - 360 000 
compound iterations (5 x overrelaxation + 1 x heat bath). All of the simulations 
in Table [I] are multicanonical, i.e. the probability distribution has been modified to 
enhance the tunnelling between the broken and symmetric phases. For technical de- 

and to ref. HIP 



tails, we refer to Sec. 



which includes a detailed description of the 
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Figure 7: The "equal weight" probability distributions of (H\H\-\- B.\B.2)jT 2 (without 
the subtractions in Eq. ( |B.7| )), measured from (3 W = 4/(g^aT) = 12 lattices. 

application of the multicanonical method to the MSSM (albeit with only one SU(2) 
Higgs doublet). 

At this point in parameter space the transition is relatively strong, as can be seen 
from the probability distributions of HlHi + H2H2 in Fig. [7]. The distributions mea- 
sured have here been reweighted to a temperature which gives equal probabilistic 
weights to the symmetric and broken phases ( "equal weight histograms" ) . 

5.4. Results for various physical observables 

Critical temperature: Tuning the temperature so that the "order parameter" prob- 
ability distributions have equal weights in symmetric and broken phases (see Fig. [F]) 
gives a good definition for the transition temperature at finite volumes. We use this 
definition in what follows. Other definitions can be found in |[L9|| ; at the infinite volume 
limit all of these give identical results within statistical errors. 

In Fig. |S| the equal weight temperatures are shown for all lattices in Table [TJ. For 
each lattice spacing we perform an infinite volume extrapolation linear in 1/ (volume). 
As seen in Fig. J8|, the slopes of the fits appear to behave somewhat unsystematically 
This is caused by the statistical errors of the individual measurements: the differences 
of the transition temperatures at different volumes are small enough that a constant 
(in 1/V) fit at each of the lattice spacings could have been used. 

On the other hand, we note a clear lattice spacing dependence between (3 W = 8 and 
the other sets of data. Obviously here the lattice spacing a = 4:/((3 w g^T) ~ 1.2/T is 
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Figure 8: The transition temperatures at different j3 w 's and volumes, plotted against 
inverse volume. The transition temperatures have been determined with the "equal 
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so large that the subleading corrections start to be significant. Thus, we use only the 
infinite volume extrapolations at (3 W = 12, 16 and 24 to estimate the continuum limit 
value, with the result 

T c = 84.3 ± 0.3 GeV. (5.11) 

This can be compared with the 2-loop perturbative value 88.5 GeV (see Figs. [|, ||). 
Indeed, the difference between the perturbative and the non-perturbative results is 
much larger than the difference between the results from different lattice spacings. 
This behaviour agrees qualitatively with our previous study ||19|. 



Triple point: In the perturbative phase diagram in Fig. |j, there is a "triple point" at 
fhu m 74.3 GeV, T c « 83.4 GeV, with v/T ~ 1.24. We have determined the triple point 
location also with lattice simulations at small volumes, using f3 w = 12, volume 12 3 (see 
Fig. |9p, and (3 W = 16, volume 16 3 . It is technically difficult and very time-consuming to 
perform multicanonical simulations at the triple point, and we did not attempt to do an 
infinite volume extrapolation here. However, the results obtained from the two lattices 
agree reasonably well with each other, and we can make an estimate for the triple point 
location, T triple « 77 ± 1 GeV, m^tripie ~ 75.3 ± 0.5 GeV, with an expectation value 
(see below for its determination) v/T ~ 1.65 ± 0.10. We emphasize that these values 
are just rough estimates; the lack of an infinite volume extrapolation can be significant 
here, since the SU(3) [/-field is very light at this point. Nevertheless, a comparison 
with the perturbative values is shown in Fig. |5|. The deviation from the perturbative 
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Figure 9: The joint probability distribution of H\Hi + H2H2 and UW at the triple 
point, measured from a /3 W — 12, 12 3 lattice. The transition symmetric^broken U 
is weak, as can be seen from the small suppression between the peaks. This means 
that it is relatively easy to end up in the phase with broken U, if one goes too close 
to the triple point. On the other hand, there is a suppression of ~ e -20 between the 
symmetric and broken H peaks, signalling a strong transition. 

triple point matches the behaviour at % = 65 GeV, where we have a much better 
control of the systematics. 

Order parameter discontinuities: Going back to rhu = 65 GeV, the discontinuity 
Av/T = [2J2i AHjHi/T 2 } 1 / 2 can be quite precisely measured from the probability 
distributions (Fig. [j]), by determining the positions of the symmetric and broken phase 
peaks separately. The results for all the volumes are shown in the left panel of Fig. [H| 
We see that for each lattice spacing, v/T remains practically constant over the whole 
range of volumes. Nevertheless, we again make an extrapolation linear in 1/V to the 
infinite volume limit. 

The lattice spacing dependence is relatively pronounced, and the discontinuity Av 
becomes smaller as the continuum limit is approached. In the right panel of Fig. [LT] 
we show the infinite volume extrapolations of v/T against the lattice spacing. The 
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Figure 10: Left: the finite volume dependence of (v/T)^, together with infinite volume 
extrapolations. Right: the continuum extrapolation of (v/T) c . Note that in units 
of T _1 , we are able to operate on relatively large lattice spacings a, since the most 
ultraviolet fluctuations of the original 4d theory have been removed by dimensional 
reduction. 



data is extrapolated to the continuum limit using a polynomial fit of the form Co + 
C2(aT) 2 + C3(aT) 3 . A priori, there is no reason why a term linear in a should not be 
there; however, when it is included in the fits, v/T invariably becomes larger as a is 
decreased at small a. While in principle possible, this kind of a behaviour does not 
seem very plausible, especially in view of the fast apparent convergence of the data as a 
is decreased in Fig. [TO. Thus, we effectively force the term linear in aT to be > in the 
extrapolation. This behaviour persists if we drop the largest lattice spacing aT =1.2 
from the analysis; in this case we obtain a good fit with the ansatz Co + C2(aT) 2 , with 
almost the same value and errors for c$. Thus, in the continuum limit, we quote our 
result as 

(v/T) c = 1.02 ± 0.05 , (5.12) 

where the error includes both the statistical errors and the ambiguity due to different 

extrapolations. This can be compared with the 1-loop perturbative result (v/T) c « 

0.47, and the 2-loop perturbative result (v/T) c ~ 0.85. See Figs. [], [TO] for comparisons. 

In addition to the discontinuities at the transition point, we have also determined 



the Higgs condensates at temperatures close to T c . In the first panel of Fig. [12| we show 
how the bilinears H 2 H 2 , R = KeHlH 2 , I = Im H1H2 behave as we go through the 
transition. Since in our simulations tan/5 = 12 is large, HlHi is practically constant 
on the scale of the plot. We note that / remains almost constant, too, and has at most 
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only a very small jump at the transition; it is non-zero overall because of the small 
imaginary parts in the effective parameters, particularly m\ 2 (T), in Eqs. (|5.4|) - (|5.10 ). 



Latent heat: The latent heat of the transition is closely related to the order pa- 
rameter discontinuities of the 3d theory. In general, the latent heat is determined 

where Af = /symmetric - /broken is the free energy density difference between the symmet- 
ric and broken phases, and the derivative is to be taken along the metastable branches. 
More concretely, L measures the discontinuity in the energy density. 

In our 3d theory with the parametrization in Eqs. (|5.4j )- (|5.10|) , explicit temperature 



dependence appears only in the Higgs mass parameters. Thus, following |53| , |55], the 
latent heat becomes 

l = rW^A^m +£ hIh*^ + (h\h 2 ± mm + H,.)\ , 

c \ dT T 2 fr[ l dT T 2 V dT T 2 //' 

(5.14) 

Wnere ZA y ■ / yj broken \ " / symmetric • 

We observe that the latent heat behaves numerically very much like the discontinuity 
in v/T, but with slightly larger statistical errors, since more condensates are involved. 
We do not show a separate figure. It should be noted that the WU-term in Eq. (|5.14fl 



is significant, despite the fact that the £7-field remains "unbroken" at all temperatures 
at this rhjj- This is because (UW) is somewhat smaller in the phase where the SU(2) 
Higgs fields H \ are "broken" . 

We use similar infinite volume and continuum extrapolations for the latent heat as 



for the condensate v/T in Fig. |10[ The final result (at our specific parameter point) 
becomes 

L/T c 4 = 0.42 ± 0.03 , (5.15) 

which is substantially larger than the perturbative 2-loop value L/T^ w 0.26. The 



difference is quite large, but it is again qualitatively similar to the one observed in [19 . 



Surface tension: We measure the tension of the interface between the symmetric 
and broken phases (also called the phase boundary) with the histogram method. The 
surface tension a equals the additional free energy/area carried by the interface. This 
extra free energy suppresses mixed phase configurations by a factor ex exp(— aA/T), 
where A is the area of the interfaces. This causes a characteristic "valley" between the 
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Figure 11: Left: the finite volume dependence of the surface tension. Right: a contin- 
uum estimate. We should stress that our determination of cr/T c 3 is not meant to be 
nearly as precise as the determination of (v/T) c (see the text). 

symmetric and broken phase peaks in probability distributions of order parameter like 

quantities; see, for example, Fig. |7|. 

In the histogram method one uses the mixed phase suppression to measure the surface 

tension. Assuming that the interfaces are perpendicular to the X3-axis direction, a can 

be obtained from 

a 

T 



1 



In 



2L1L2 P u 



as V — > 00 



(5.16) 



Here P, 



max. mm 



are the maximum and minimum of the peaks of the probability distri- 
bution, and Li = N^a are lattice extensions in physical units. The factor 2 appears 
in front of L1L2, because there are two interfaces in a box with periodic boundary 
conditions. 

In the actual analysis we measure a from histograms of H\Hi + h T \h T 2- These are 
reweighted to a temperature where the peak heights are equal, which simplifies the 
analysis. We also use a finite volume scaling ansatz similar to our earlier work |19[ |, 



a 
j>3 



2(LTf 



p 

1 r 1 



In — — - H — In LT + const. 

■L min ^ 



where we have employed the fact that our lattices are all cubic, L\ = L2 = L% - 
In Fig. [IT] we show the rough infinite volume extrapolations of <r/T c 3 from j3 u 



(5.17) 



12, 



16 and 24 lattices. (For this observable, (3 W = 8 gives a surface tension larger by about 
a factor of 5, and we do not include it in the analysis.) The continuum limit estimate 
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of a becomes now 

ajTl « 0.010 ± 0.05. (5.18) 

The 2-loop perturbative value is cr/Tj? « 0.017, which is larger than the lattice value. 
However, we stress first of all that our lattice determination is quite rough here, as 
we have only used cubic volumes, not the cylindrical ones often employed for getting 
good infinite volume extrapolations (see, e.g., JBJ). Second, we should also remember 
that in the perturbative estimate no account is taken of the effects related to derivative 



terms discussed in N76L W% W& 79, 351, which would decrease that value. Nevertheless, 



the situation is certainly different from the case studied in |L9] where even the raw 



perturbative value was smaller than the lattice value. On the other hand, the situation 



is qualitatively similar to the case of the Standard Model [55], where the transition is 



of a relatively weak strength compared with |I9], as it is here too. 



Correlation lengths: Next, let us explore spatial correlation lengths around the 
transition temperature. As usual, their inverses are called "screening masses". The 
masses are measured from the 0-momentum correlation functions 

C a (x 3 ) = ^(O„(y,0)O a (x,x 3 )) ocexp[-m a x 3 ], (5.19) 

x,y 

where x, y are in the transverse (xi,X2)-plane, and the gauge invariant operator O a (x) 
is one of the following: 



C even SU(2) scalars: 


s:( X ) 


= Hl(x)H a {x), R(x) = ReH\(x)H 2 (x), 


C odd SU(2) scalar: 


I(x) = 


Im Hl(x)H 2 (x), 


SU(2) vectors: 


V a %x) 


= 1mHt{x)U?(x)H a {x + i), 


SU(2) 0++ glueballs: 


G w (x) 


= l-±TrPf 2 (x), 


SU(3) scalar: 


S s {x) -- 


= W(x)U{x), 


SU(3) vector: 


Vf(x) - 


= lmW(x)U l s (x)U(x + t) 1 


SU(3) 0++ glueballs: 


G s (x) -- 


= 1- |ReTrPf 2 (x). 



Here i — 1, 2; a — 1, 2 labels the two SU(2) Higgs fields; Uf and Uf are the standard 
lattice SU(2) and SU(3) gauge links (denoted by U% it U^ in Eq. Q); and P&, Pf 2 are 
(x±, x 2 )-plane plaquettes constructed from these links. In order to reduce statistical 
noise, we use recursive blocking and smearing of the gauge and Higgs fields along 
(xi,X2)-planes, and construct operators and measure all of the correlation functions 
from the blocked fields at each blocking level. The blocking is repeated up to 4 times. 
For details of the recursive blocking procedure we employ here, we refer to [|l9j (see 
also H|). 
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Figure 12: Top left: Higgs bilinears, subtracted according to Eq. ( |B.7| ), across the 
transition at (3 W = 16, volume = 32 3 . The vertical dotted line shows the location of 
T c . Top right: light SU(2) scalar and vector inverse correlation lengths, m/T. Bottom 
left: heavy SU(2) scalar, SU(2) glueball and SU(3) scalar inverse correlation lengths. 
Bottom right: inverse correlation lengths related to the operators R, I. All the numbers 
in the bottom panels have relatively large systematic errors (see the text), and could 
almost be compatible with each other. 

The measurement of the correlation lengths in an interacting theory is complicated 
by the fact that all the operators in a given quantum number channel in general couple 
to the same set of physical states. For example, we can expect that all of the scalar 
operators above, including the ++ glueballs, will in the limit 23 — > cxd yield the 
screening mass of the lightest scalar state. On the other hand, in the real world 
the correlations may behave at intermediate distances in a different way, and to fully 
resolve this behaviour one usually measures the cross-correlation matrix of a large set 
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of operators in a given quantum number channel, and diagonalises it. 

However, in our case when only rough qualitative accuracy is needed, it turns out that 
this is not necessary: since we use a large tan/3 = 12, the active Higgs component in the 
transition projects almost completely to H 2 , and correspondingly the light scalar and 
vector states couple strongly only to the operators Sf and V 2 ™, respectively Moreover, 
the "coupling" (or overlap) of the SU(2) pure gauge and the whole SU(3) sector to 
the SU(2) scalar Higgs sector is weak. All in all, this implies that we can obtain the 
heavier scalar and vector screening masses just by measuring the exponential fall-offs 
of the corresponding correlation functions at intermediate distances; the accuracy of 
this approach is quite sufficient for our conclusions. 

On the second panel of Fig. [T2] we show the light SU(2) scalar and vector masses, 
and on the third and fourth panel the heavier masses. The lightest mass scale at the 
transition is an order of magnitude smaller than the heavy scales. This makes the 
condition in Eq. (|3.2j ) very difficult to meet in practice, and in our case the first part 
a <C £ min is barely true. This circumstance also makes the measurement of the heavy 
masses difficult, since the correlation functions vanish into noise after a few lattice units, 
which explains the large errors in the data. However, for our purposes the accuracy 
obtained is sufficient. 

We may now note first of all that while the large mass scales make the extrapolation 
to the continuum limit delicate (necessitating several lattice spacings, as we have), 
experience from QCD |HJ indicates that the mass scales < 3T are still far from being 
too large for us to have any concerns about the applicability of dimensional reduction, 
used in the construction of the 3d effective field theory. The integration out of the 
Matsubara zero modes of the temporal components of the gauge fields is clearly more 



critical, but as we have discussed in |49| , we expect even that to be reasonably under 
control. In any case, those masses are still much above the lightest ones in the sys- 
tem (<^C T), which determine the non-perturbative thermodynamical properties of the 
transition. 

Furthermore, the screening mass spectra measured add further confidence to two 
statements we have made on other grounds before: (1) The correlation length related 
to / is short, and thus not at all "critical". This means that we are far from the 
possibility of spontaneous C violation. (2) The heavy SU(2) scalar correlation length 
is much shorter than the light one, and thus again far from "critical" . This means that 
it is only the "light" combination of the two Higgs doublets which is really "dynamical" 
at the transition point, and any effects related to the other one are suppressed. 
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6. The properties of the physical phase boundary 

We now turn to the study of the properties of the phase boundary. We have outlined 
the measurements to be carried out, as well as some caveats in them, in |49| . 



In order to study a phase boundary, we have to make sure that there really is 
one on the lattice.^ This can be achieved by restricting, say, the volume average of 
h 2 = H\Hi + H2H2 to a narrow band around (/i 2 ymmetr i C + ^broken)/2- Provided that 
the volume of the system is large enough, this guarantees that the system will always 
remain in a broken + symmetric mixed state, with corresponding phase boundaries, or 
interfaces. 

In the case at hand the interfaces are rather thick, and because of the periodic 
boundary conditions, there will be two interfaces spanning the lattice. This makes it 
advantageous to use cylindrical lattices: because of the surface tension, the interfaces 
will be preferentially oriented along the smallest cross-sectional area across the lattice, 
making them well separated along the longest lattice direction (23, say). This has the 
further advantage that we always know the orientation of the interfaces. Our interface 
simulations were made on j3 w = 8, 12 2 x 96, and f3 w = 16, 24 2 x 192 lattices, using up 
to 450 000 compound update sweeps per lattice. 

6.1. Observables 

We study the interface properties by looking at the Higgs bilinear operators in Eq. Q3.3|) 
as follows: first, for each configuration, we average the bilinears across the (21,22)- 
plane, so that we obtain them as functions of the 23-coordinate. Then we average 
over all configurations in two ways: (1) the bilinears are measured as functions of the 
distance from a certain reference point, and (2) one bilinear is measured as a function 
of another. Let us look at these cases separately. 

(1) Interface profiles: The Monte Carlo simulation method described does not 
specify the location of the interfaces along the 23-direction. However, in order to 
measure the profiles of various observables across the interface, we have to find a 
reference point in the 23-direction, configuration by configuration (i.e., we have to 
remove the zero mode). Some care has to be taken here: for example, one may locate 
the 23-value where h 2 , averaged over x\ and 22, reaches the value half-way between the 

1 In principle, standard (multicanonical) simulations at the transition temperature could be used, 
since phase boundaries appear there in the "tunnelling configurations" containing regions of both 
phases. However, a lot of effort is wasted, since phase boundaries exist only in a small subset of the 
total of all configurations. 
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symmetric and broken phase ones. However, since h 2 (as any other local quantity) has 
large fluctuations, this kind of a sharply defined location will cause the configuration 
to be shifted such that the natural fluctuations across the interface may be summed 
"in phase" , which tends to distort the profile. The problem can be avoided by using a 
"softer" filter function. In this work employ the fact that there are two interfaces on a 
periodic lattice: thus we can find the "symmetry point" 23,0 of the profile by taking a 
Fourier transform (N 3 is the extent of the lattice in the direction of x 3 ), 

C = J2h 2 (x 3 )e l2nX3/N3 = Ae l27TX3 ' o/N \ (6.1) 

Thus £3,0 = N 3 /(2tt) x (polar angle of C) . Configurations are then superimposed 
around this point. 

In Fig. [13] we plot the profiles of the bilinears H1H1, H2H2, R = Reifji^ and 
/ = ImHlH2, measured from the (3 W = 16 lattice. (We observe no qualitative change 
compared with the /3 W — 8 lattice.) Only half of the profiles (one interface) are shown. 
Besides the obvious one in magnitudes, it is difficult to see qualitative differences 
between H\H^ H2H2 and R across the interface. 

As for the C odd condensate /, its magnitude is even smaller and the errorbars cor- 
respondingly larger. For clarity, we have smoothed I(x 3 ) in Fig. [IB] by an approximate 
Gaussian smearing: 

I{x 3 ) <- i [I{x 3 -a) + I(x 3 ) + I(x 3 + a)} , (6.2) 

and repeated 4 times. Without this additional smearing it would be difficult to see any 
structure in the plot. After the smearing, we observe that I increases slightly when the 
broken phase is entered. The overall negative value of I is due to the small imaginary 
part of the m 2 2 (T) parameter, Eq. (|5.6|) (see Appendix |A.1|) . 



From Fig. |13] we see that the interface is rather thick: if we fit H\Hi to a function 
of the form a + 6tanh((x3 — c)/L), we obtain L ~ 9/T. However, one should bear in 
mind that in 3 dimensions the apparent interface thickness suffers from a logarithmic 
divergence as the area is increased (see, e.g., |82]|). A natural way to resolve this 



arbitrariness is to consider the interfaces on physically relevant length scales. For our 
case, the longest correlation length £ max at the transition is ~ 6/T (see Fig. O), so 



that the interfaces in Fig. 13 correspond to a cross- sectional area ~ (2£ max ) 



2 



(2) H1H1 vs. H2H2 and I vs. R across the interface: Instead of plotting the 
bilinears as functions of x 3 , it can be more illustrative to consider the behaviour of 
one condensate as a function of another. This way there is no logarithmic divergence 
visible, and no need to locate the interface on the lattice. 
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Figure 13: The profiles of the operators in Eq. ( |3.3|) as a function of the spatial coor- 
dinate across the wall. The cross-sectional area here is ~ (14/T) 2 . The subtraction in 
Eq. ( |B.7| ) has not been carried out. 



In Fig. fl we show H 2 H2 as a function of H\Hi, and I as a function of R. In the 
former case, there is a small but statistically clear deviation from the straight line 
between the symmetric and broken phases. We show the same data in Fig. |T5| using 
the definition in Eq. (|3.5|). 

We observe that in terms of t>j = y2AHJHi, the deviation from a straight line is 
0(a. few x 10 -2 ). We also observe that C violation has some structure: the operator / 
tends to saturate close to the broken phase value considerably faster than R (or v/T). 
However, we do not observe any significant amplification of / inside the interface, and 
thus no sign of transitional C violation. 

The orders of magnitude observed for the variation in / and tan/3 agree roughly 
with perturbation theory |)7], ^, |3(| and, in the case of tan j3, with previous 4d simu- 
lations [Q, although in most cases the parameter values were not identical to ours. 
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Figure 15: Left: behaviour of the effective tan /3, Eq. ( |3.5| ), across the interface. We note 
that the "dynamical" value differs only a little from the zero temperature input value 
tan/? ps 12.0. Nevertheless, the effective tan/? is not constant across the interface, but 
varies by a relative amount < 5%. Right: the same data, plotted as (v /v broken) (tan (3 — 

tan Awoken)- 



34 



7. Summary and Conclusions 

In this paper, we have studied the electroweak phase transition in the MSSM, with 
particular attention on CP violation in the background ("vacuum") configuration, as 
well as on the strength of the phase transition. 

The method we have used is based on 3d effective field theories, and their non- 
perturbative study. At finite temperatures around the electroweak phase transition, 
the thermodynamics of the MSSM can be represented by a theory containing two SU(2) 
Higgs doublets and one SU(3) stop triplet. Despite its complexity, we have demon- 
strated that this theory can be studied in a controlled way with lattice simulations. 

The phase diagram of this theory is non-trivial, involving a phase where CP is (even) 
spontaneously violated, as well as a phase where the U(l) symmetry corresponding to 
a massless photon is broken. We have studied the phase with spontaneous CP violation 
in some detail. We have found that for the parameter values allowed by the MSSM, 
one does not end up in this phase close to the electroweak phase transition. In more 
general two Higgs doublet theories this could happen, with potential implications for 
baryogenesis. 

We have then studied the electroweak phase transition at physical parameter values, 
in particular m# ~ 105 GeV, not ruled out for the MSSM. We observe very clearly the 
feature familiar from our previous MSSM study [|19| that the transition is significantly 



stronger than in 1-loop perturbation theory, and even stronger than at 2-loop level, due 
to the fact that the critical temperature T c is lower. Let us note that the situation is 



different from that studied with 4d simulations in ||62|| , where the transition was quite 
strong (rriH « 45 GeV), and good agreement even with 1-loop perturbation theory was 
found. We do not consider there to be any discrepancy, however, since all our previous 
experience is that perturbation theory works the better, the smaller the Higgs mass. 

At the point of our present study, we observe (Fig. [10]) that the transition is strong 
enough for baryogenesis, since v/T ~ 1.0. (Based on analytic estimates, one would 
expect that v/T has to be somewhat larger than 1.0 []55| , but a dynamical lattice com- 
putation suggests that 1.0 should be enough |Q). However, this concerns a particular 
parameter point, and it is important to ask how much room there is around it. 

To get a comprehensive estimate, let us now go to the point of the strongest possible 
transition, i.e. the triple point shown in Fig. |5|. We then vary various parameters and 
use 2-loop perturbation theory to compute v/T. The results are shown in Fig. |I| We 
find the rather remarkable fact that the results are almost independent of the Higgs 
sector parameters rattan/?. This is so because at the triple point the properties of 
the transition are dictated by the stop sector. 
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Figure 16: The Higgs mass m# and broken phase expectation value v/T at the triple 
point (see Fig. |]) according to 2-loop perturbation theory, as a function of various 
quantities. The vertical lines indicate the location of our reference parameter values. 
The Higgs mass is computed with the 1-loop relations specified in |49| (corresponding 
to f84|). Note that the combination relevant for large rriA, A t — fi* cot j3, is in most 
cases close to A t , since we keep /i = i200 GeV fixed and cot/? = l/tan/5 is small. In 
the values of ran we expect an uncertainty of a few GeV, due to unimplemented 2-loop 
zero temperature Higgs mass corrections, 2-loop dimensional reduction corrections, and 
perhaps also the fact that explicit CP violation has not been rigorously treated in the 
zero temperature relations we employ here ||9) (for recent discussions, see [Q 
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On the other hand, the transition weakens rapidly away from the triple point (see 
Fig. |5]). Our lattice results provide a strengthening effect which can partly compensate 
for this. Nevertheless, one needs to remain close to the triple point in any case, for 
instance rhy ~ 65... 77 GeV for the parameters employed in Fig. |. The perturbative 
range would have been rhjj ~ 69... 74 GeV. 

Apart from v/T, we have also measured other important characteristics of the phase 
transition. The values of the latent heat L/T* and surface tension <t/T c 3 allow us to 



discuss the real time history of the phase transition. Estimates such as in [§5, |86], 87, 
[I7J, [H| lead to the conclusion that the latent heat is probably large enough to reheat the 
system back to T c after the bubble nucleation period |35|, since L/T* > 8(<r/T c 3 ) « |^7J . 



(In fact, the system looks very much like the case B studied in p7[ , but the physical 
friction is orders of magnitude larger |33j than assumed in |T7| based on the literature 
available at the time, and the physical velocities are therefore smaller than in [J87 , 
< 0.1 |33], |34]|.) The small bubble wall velocities before and particularly after reheating 
(when they are ~ 0.001) may lead to enhanced baryon number production according 



to the standard computations [24, 25, 26 



Finally we have studied the properties of the phase boundary, or bubble wall, at 
the physical transition point. We have determined the profiles corresponding to tan/? 
and to the C violating phase angle cos0 numerically, and excluded spontaneous (also 
called transitional) CP violation within the phase boundary, too. Explicit CP violating 
effects in the Higgs background are non- vanishing but small, even if the explicit phases 
are of order unity, because they are suppressed by effective couplings of the type g^T 
over the heavy mass scale >(m^ + 0.5T 2 )2. The profiles we have determined could 
in principle be used as the semiclassical background entering the actual baryogenesis 
computations [p4f-||31~ll,f37 



In summary, from the point of view of the non-equilibrium constraint, there is 
some parameter space available for electroweak baryogenesis in the MSSM. Our non- 
perturbative results agree well with the ones in p0[, based on 2-loop perturbation 



theory and our previous lattice results [11|, and allow for a strong transition even 
for a Standard Model like Higgs mass m# ~ 115... 120 GeV if tuq > a few TeV and 
A t /rriQ > 0.5 (see Fig. [16]). On the other hand, we find also a small value of rriA < 120 
GeV to be acceptable close enough to the triple point, even though away from it large 
values are favoured. Small values of m^ make the MSSM look less like the Standard 
Model and relax the experimental constraint on m H fLOf , |63|]-|)(J, allowing perhaps 
for somewhat smaller itlq as well. Thus electroweak baryogenesis continues to be a 
viable scenario, besides for instance those based on Majorana type neutrino masses, if 
at the same time quite strongly constrained. 
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Appendix A. Integrating out the heavy Higgs direction 

We review in this Appendix how the effective theory in Eq. ( |2.1| ) can be simplified by 
integrating out a linear combination of the Higgs doublets, if we are not interested in 
C violation but only in the strength of the phase transition. We have discussed the 
procedure previously in Sees. 6,7 of [|l(| and in Sec. 3.1 of J49| . We complete those 



results here by allowing for complex parameters (explicit CP violation), as well as by 
having a light dynamical stop. We work at 1-loop level. 

It should be noted that contrary to the case at zero temperature, integrating out 
a linear combination of the Higgs doublets is reliable even for small values of m^, 
because thermal corrections increase the effective mass of the degree of freedom that 
is integrated out (see below). 

A.l. Phase redefinition 

The starting point is the effective theory in Eq. ( |2.1| ). We take first a trivial step, 
removing one extra phase from the parameters in order to simplify the notation. Indeed, 
if m1 2 (T) = \m\ 2 {T)\ exp(i0 12 ), then we can make a field redefinition 

^^H^ 12 , H 2 ^H 2 . (A.l) 

As a result, the real parameters in Eq. (|2.1|) remain unchanged, but the five complex 
parameters change as 

4 M) (T) = \m\ 2 {T)\, (A.2) 

(A3) 
(A.4) 
(A.5) 
(A.6) 
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(new) 
7l2 


= 7i 2 e- i012 , 


^ (new) 
A 5 


= A 5 e- 2 ^ 12 , 


\ (new) 

A 6 


= A 6 e"^ 2 , 


•i (new) 

A 7 


= A 7 e"^ 12 . 



We leave out the superscripts "(new)" in the following, with the understanding that 
after each step, the new parameters are denoted with the same symbols as the old ones 
before it. 

Note that if 712 , A5, Xq, A7 are small and there is no spontaneous CP violation, 
then this field redefinition directly determines the phase angle of (H\H 2 ). For in- 
stance, for the parameter value in Eq. (|5.6| ) with T < 100 GeV, we get m\ 2 {T) ~ 
-1200 e iom7r GeV 2 , and consequently H\H 2 ~ \HlH 2 \ e - iomw , and \mH\H 2 < 0; see 



Fig. 14 



A. 2. Diagonalising the mass matrix 

Next we want to define new fields as linear combinations of H\,H 2 , such that the term 
mixing the two directions, ~ m\ 2 (T)H\H 2l vanishes at tree-level (1-loop corrections 



can still induce a mixing and this effect shows up below). Following 

H 1 = cosa h + sina H, 



H 2 



sina h + cosa H. 



The angle a is chosen so that 
2m 2 2 (T) 



tan2a 



m 2 {T)-m 2 {T) 



sin2a 



2m 2 2 (T) 



m\{T) - m|(T)) 2 + 4mf 2 (T) 



we write 

(A.7) 
(A.8) 

(A.9) 



It should be noted that in the practical case considered in Sec. |5|, the large value of 
tan/? implies that a ~ 7r/2, which means that the light field h is almost in the direction 
of the original H 2 . 

After the rotation, the quadratic part of the scalar potential is 



V, 



m T 



'2 — "<>U 

where the new mass parameters are 

1 



(T)U^U + m 2 h {T)h)h + m 2 H (T)H^H, 



m 2 h (T) 

m 2 H (T) 



m\{T) + m 2 2 (T) - J{m\{T) - m 2 {T)f + 4mf 2 (T) 



mi(T) + m 2 2 (T) + J(m 2 (T) - m 2 {T)Y + 4mj 2 (T) 



(A.10) 

(A.ll) 
(A.12) 



The stop mass parameter m^T) does not change from the value in Eq. (|2.1|) . 

The scalar couplings are modified as follows. The stop self-coupling Xu does not 
change. Denoting the quartic scalar potential related to the Higgses by 



V/Lffi 



ggs 



+ 



jxrfUtfh + l2 U ] UH^H + [ ll2 U ] Uh)H + H.c. 
Ai(/A) 2 + X 2 {H^H) 2 + X 3 rfhH^H + XjJHrfh 
X 5 {h ] Hf + X 6 rfhtfH + X 7 H^HrfH + H.c. 



(A.13) 
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we get 



(new 

7i 

(new 
72 

t-, (new 

Re7} 2 

T (new 

I1II7I2 

ImA 7 ncw 



71 cos 2 a + 7 2 sin 2 a — Re 712 sin 2a, 
71 sin 2 a + 72 cos 2 a + Re 712 sin 2a, 

7: (71 ~ 72) sin 2a + Re 712 cos 2a, 

Ini7i2, 

Im A5 cos 2a + - Im(A6 — A7) sin 2a, 

— Im A5 sin 2a + Im A6 cos 2 a + Im A7 sin 2 a, 
Im A5 sin 2a + Im X 6 sin 2 a + Im A7 cos 2 a. 



(A. 14) 
(A.15) 

(A.16) 

(A.17) 

(A.18) 

(A.19) 
(A.20) 



The Higgs self-couplings A1...A4, together with the real parts Re A 5 ... Re A 7 , on the other 
hand, are related by the matrix in Eq. (6.21) of |l6f . 



A. 3. Integrating out the heavy direction 



In Eq. QA.9Q the angle a has been chosen such that the field h is light, as can be seen 
from Eq. (|A.11|) . Then the heavy field H can be integrated out. Indeed, the expansion 
parameters related to this integration are 



9 2 W T 



XT 



4irm H (T)' Anm H (T)' 



(A.21) 



which are small close to the phase transition. This is because one of the eigenvalues 
of the Higgs mass matrix must be very light at the point of the phase transition, 
m1(T) ~ (g^T) 2 , so that the other eigenvalue is equal to the trace of the mass matrix, 
given in Eq. 



m 2 H {T) ~ m\{T) + m 2 2 (T) >m 2 A + 0.5T 2 . 



(A.22) 



When H is removed, the resulting theory is just the same as studied in [17], [L9 

£ 3d = ^Tr (^ j + {DtU)\D' i U)+rr^(T)U^U + X u (lPU) 2 

+ \TiF 2 j + {D^h)\D^h)+m 2 h {T)h)h + \ h {h)h) 2 + 1 U^Uh)h. (A.23) 
2 J 



At 1-loop level, the new couplings are: 

9lT 



2(new) 



9 2 J'i 



A87rm H (T) 



(A.24) 
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9 2 s {new) = 9l (A.25) 

m 2(new) (T) = m 2 (T) _ _L (2Ag + Xi )m H (T)T, (A.26) 

m 2 J now) (T) = ml{T) - ±- l2 m H (T)T, (A.27) 

Ai new) = Ai - - T 7^ (X 2 3 + A3A4 + ^A 2 + 2|A 5 | 2 + 12Re(A 6 - A 7 )Ag),(A.28) 

7 (~w) = _ r ( 2A , i72 + A + 4 | |2 + 12 Re(A<i _ At) . ■) (A30) 

A. 4. 2-loop mass parameters 

Finally, let us recall how the results above would change by a 2-loop integration out 
of H. From the practical point of view, the most important effects are in the mass 
parameters f52f . After the integration, the renormalized mass parameters in the MS 
scheme can be written as 



m 



>(T) = m 2 (T)+ / 4 +9A ^2_ m 2_3 7 2 + 8 ^ ln ^ (A>31) 



,, (T) = ml(T) + 7 ^(8^ 4 + ^A^ s 2 -16A^-2 7 2 + 3g 2 wl ) In ^, (A.32) 



2(new) 

m 



(An) 2 v 3 H 



where m1(T),mjj(T) are the 1-loop results in Eqs. ( |A.26| ), ( |A.27|) . Thus a 2-loop 



computation amounts to a determination of the expressions for A^, Au ~ a few x 



T |52|, O, pTL [731; see Sec. IO for a discussion of the status of such computations. 



Appendix B. Lattice counterterms 

We collect here the lattice counterterms needed in Sec. 13.11. The derivation of the 



counterterms proceeds as in p31 BR 1811 , and a major part of the results can be extracted 



from there. However, some new parts are needed too, because there are now two SU(2) 
Higgs doublets in contrast to just one. 

The most non-trivial 2-loop changes can be obtained as follows. In the contributions 
proportional to g 4 , we have to replace T by Y^iTi in Eq. (E.4) of f89j , where i runs 
over all the fields (fundamental or adjoint) interacting with the SU(iV) gauge fields, 
and Ti — I'm the former case, N in the latter. In the present case of two fundamental 
doublets, one thus simply needs to put J2i T% — > 2 for the SU(2) case g = g w . To obtain 
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the g 2 X, g 2 7-terms, we replace m 2 d by the trace of the scalar mass matrix, computed in 
the appropriate Higgs background, in Eq. (E.5) of pSfl . Finally, the numerical factors 
in the terms of types A 2 , 7 2 have to be computed by hand. 

The bare parameters appearing in the lattice action are then of the form 



m 



LB 



m 2 (T) + Sm 2 . 



(B.l) 



where m 2 (T) are the MS scheme parameters at a scale /i. The results for the counter- 
terms 5m 2 are: 



dm 2 ; 



■^(J^ + 8A c/ + 2 7i + 2 72 )t 

rp2 



16vr 2 



64 

8 at + ~j Xu9 * + 3 ^ 71 + 72 ) 

6 
aj2 



-16A 2 ; - 2( 7l 2 + 7 | + 2| 7l2 | 2 )) An — + 0.08849") 



+19.633^ s 4 + 12.362A^ 2 + 1.7384( 7l + l2 )g\ 



(B.2) 



5m\ 



a 



91 



47raV2 
T 2 r/45 



6Ai + 2A 3 + A 4 + 3 7 i T 



16tt 2 



, g^t + 2 ( 6A i + 2A 3 + X ^l + 8 W - 3( 7l 2 + I712I 



5m 2 , 



-(12A 2 + 2X 2 3 + 2A^ + 2A3A4 + 12|A 5 | 2 + 9|A 6 | 2 + 3|A 7 
+5.46500* + 0.86921 (6Ax + 2A 3 + X 4 )g 2 w + 4.6358 7l ^ 2 

(l 



6 



In — + 0.08849 



(B.3) 



4vraV2^ 
T 2 r/45 



0; + 6A 2 + 2A 3 + A 4 + 3 7 2 )T 



16tt 2 



^i + ^2 + 2A 3 + K)g 2 w + 8 l2 g 2 - 3( 72 2 + I712H 



5m 2 u 



-(12X 2 2 + 2A 2 + 2A^ + 2A3A4 + 12|A 5 | 2 + 3|A 6 | 2 + 9|A 7 | 2 )) fin — + 0.08849 
-5.46500* + 0.86921 (6A 2 + 2A 3 + X 4 )g 2 w + 4.6358 72 # 2 , (B.4) 

A 3 (A 6 + A 7 + 7l2 )T 

/g 

( 2^ 6 + X ^ 9 ™ + 8 ^ 12 ^ ~ 3 7i2(7i + 72) 



T 



16tt 



'9 

:2 { 

-3(2AiA 6 + 2A 2 A 7 + (A3 + A 4 )(A 6 + A T ) + 2A,( A;, 1 A r 
+2.6076(A 6 + X 7 )g 2 w + 4.6358 7l2 £ 



In — + 0.08849 



(B.5) 
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Here S = 3.175911535625 and a is the lattice spacing. 

The continuum operators in Eq. ( |3.3| ), on the other hand, are obtained as 



(Wu) 



J"2 



MS, a 



{frV) 



lattice 



^T^- (b - + ° M7% 



j> 2 



MS, 



(HjHj) 



,i* 



lattice 



+ 

E 



a/j, 



, +-^L fin— + 0.66796) 

2iraT (An) 2 V afi J 



In practice we choose to discuss MS parameters with \i = T, so that 



6 



6 



a\i aT 2 



S s 2 &. 



(B.6) 
• (B.7) 



(B.8) 



Appendix C. The C violating phase in perturbation theory 

We collect here the details related to the discussion outlined in Sec. W7A The starting 



point is the effective potential in Eq. (|4.2p . Note that we are free to choose (3, 9 G 
(0, f),0 G (— 7r, 7r). We ignore first A t , p, and the 1-loop effects from the SU(2) Higgs 
masses m 2 Si , and present a complete parameterization for the C violating phase in that 
case. We then discuss the effect of A t , fi ^ and m 2 Si ^ 0. 



C.l. Minimization with respect to 9,(f> 

Let us assume for the moment that v,/3, or t>i,i>2 > 0, are given. We denote 

M\ 2 = m\ 2 (T) + -\ 6 v\ + -X 7 v 2 2 , 



(C.l) 



and assume first that M\ 2 ^ 0. Minimizing Eq. ( [4.2] ) with respect to <p, we obtain that 
the region for spontaneous C violation is 



A 4 -2A 5 <0, Kttj^ > 1, 

TO 



and then 



cos 6* = 1, cos i 



M x 2 2 



2A 5 wif 2 



The region where U(l) is broken is 



A 4 -2A 5 >0, (A 4 -2A 5 )^^>1-A 5 - 2 ^ 2 



121 



(C.2) 



(C.3) 



(C.4) 



121 
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and then 



2lM 2 I 

cosfl= , ' , 12 ' , I cos 01 = 1. (C.5) 

(2A 5 + A 4 )W ' ' K ' 



For A4 — 2A5 = and As^^r > 1, cos 9 and cos cj) are undetermined but 

M 2 
cos 9 cos = r J *—. (C.6) 

Elsewhere, cos# = | cos0| = 1. The special case M 2 2 = can be treated as a limit of 
these formulas. 

We will in the following concentrate on the case in Eq. ( |C.2| ). Then, the value of the 
effective potential at the minimum of Eq. ( |C3| ) is 

^0i^2)| Eq . ( C3) = 2 m ^ T ^i + 2 m 2( T )^2 + ^M + 4A2V2 + ^v\vl 
+ ±(A 4 - 2A 5 )^ 2 2 - ^(rn 2 12 (T) + ^\ G v 2 x + ^A 7 t; 2 ) 2 

T :9 3 M + v 2 2 )l - ^K(T) + \hlvlf. (C.7) 



167r *«v 1 ^ 2 7tV u " ' 2 



C.2. Boundedness 



Next, we discuss which values of the couplings naively leading to spontaneous C viola- 
tion are actually allowed from the point of view of the consistency of the theory. Let 
us first of all recall that according to Eq. (|C.2f) , 

A 5 >0, A 4 -2A 5 <0. (C.8) 

Furthermore, for the theory to be bounded from below, we must clearly also require that 
Ai, A 2 > in Eq. (|4.2|). However, this is not enough. It turns out that the most critical 
direction in the field space is where spontaneous C violation indeed takes place (since 
this means that the 2nd order polynomial in cos0, Eq. (|4.2| ), has been successfully 
minimized). The value at the minimum is given by Eq. (|C.7|) . We observe that the 



contribution in Eq. (|C.7|) effectively normalizes the values of A1...A3 in Eq. (|4.2|) . It is 
then easy to see that boundedness requires that in addition to Eq. (|C.8|) , one has to 
satisfy 



5- 



A^ < 4AxA 5 , X z 7 < 4A 2 A 

A3 + Al _ 2A5 _ ^ > _ 2] l {Xl _A ){ ^3.). ( C9) 

These will be replaced by stronger constraints below when we restrict ourselves to 
finding a C violating minimum at some finite values of v 1 , v 2 , but are nevertheless 
useful as simple relations involving the quartic couplings only. 
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C.3. Stationary point with respect to t>i,t>2 



Next, we should minimize the effective potential with respect to vi,v% in addition to 
9,(f) as has been done before, in order to express t>i,t>2 in terms of the parameters of 
the theory. It turns out that it is convenient to turn around the question: we will use 
v i, t>2 to parameterise different theories leading to spontaneous C violation, and express 
m\(T) 1 m\(T) 1 m% 2 (T) in terms of these. 

Since the potential has already been minimized with respect to 9,<p (c.f. Eq. ( |C.7[ )), 
it is sufficient to impose dV/dvi = 0, i = 1, 2. We then find that a stationary point at 
(fi,f 2 ) = v (cos/3, sin/3), with a C violating angle cose/), is obtained for given A1...A7 
provided that the mass parameters are 



m| 

^2 



m 



12 



^2 



- A6 cos cf) sin 2/3 + Ai cos 2 /3 + — (A3 + A4 
-A7 cos 4> sin 2/3 + A 2 sin 2 /? -\ — (A3 + A 4 • 
-i7^(M 2 + t» 2 sin 2 /3) 1/2 , 
A5 cos sin 2/3 + -A6 cos 2 /3 + - A7 sin 2 (5 



where we have denoted 
3 . 



G 



16tt 



9u 



0.018, H 



A^2n 



-h\ w 0.169, 



-2A 5 )sin 2 /3 


y2 ~ 7 1 ' 


(CIO) 


- 2A 5 ) cos 2 /3 


t; 2 « 

^2" + G ^ 






(C.ll) 


v 2 

^2' 


(C.12) 


M2 = l m ' 


>(T). 


(C.13) 



C.4. Local minimum with respect to V\,V2 



Not all of the stationary points obtained through Eqs. ( p,10| )-( |C.12j ) are local minima. 
The final stage is imposing this condition, which leads to some further restrictions 
on the parameters (and on the values of v allowed). Of course, the requirement of 
obtaining a global minimum in addition to a local one, would lead to still stronger 
restrictions, but for the present purpose it is enough to consider the local condition. 

As in the previous paragraph, after the minimization with respect to 6, <fi has been 
carried out, leading to Eq. ( |C.7D , it is enough to consider the potential as a function 
of v 1, v 2- The constraint is that the mass matrix Adij = d 2 V/dvidvj have only positive 
eigenvalues, i.e., 

det.M>0, TrM>0. (C.14) 

These conditions result in the following constraints: 



Ai 



\ 2 
4A, 



GT 

~2v 



>0, 



(C.15) 
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A A? 



GT 

~2v 



HT 



2JM 2 + v 2 



>0, 



(C.16) 



A3 + A4 — 2A5 — 



AeA? 
^A7 



G 



T 



< 



\l 



A Ai 



~2v~ 



A, 



4A, 



GT 



#T 



2JM 2 + u| 



(C.17) 



Note that these equations cannot be satisfied at arbitrarily small values of v/T, and 
thus spontaneous C violation can only take place at sufficiently large v/T. 

C.5. A complete parametrization 

In the previous paragraphs, we have obtained expressions for the mass parameters 
leading to spontaneous C violation, Eqs. ( |(J.10| )-( |U.12| ), but also a number of con- 
straints that have to be satisfied, Eqs. ( P.15|) - (|C.17| ). We can now present a complete 
parametrization for all the C violating states allowed by the potential in Eq. (|4.2p (with 
A t ,fi,m 2 Si = 0), such that the constraints are automatically taken care of. 

Suppose we want to have a local minimum where C is spontaneously violated 
(| cos0| < 1), at a given vev v/T, with a given tan/5 = t^/fi- Take arbitrary Ai, A2, 
A5 > 0, and G, H, M as defined in Eq. ( |U.13|) . Then there is a 4-parameter family of 
possibilities, parameterised by 



a 3 ,a 6 ,a 7 e (0,vr); p 4 G (0, 00) 



v 



provided that 

where k is the unique root in the range 



G 



> max I - ., 

T \2Aisin olq 



H 



K > 



G 



2A 2 sin a-j 



of the equation 



2,„2 



(M 2 /T 2 + k 1 sin i /5)(2KA 2 sin^a 7 - G) 2 = H 2 k 



The remaining couplings have to be chosen as 



T 



A 3 = 2p 4 + G- + 2 JAxA 



cos a 6 cos a 7 



(C.18) 
(C.19) 

(C.20) 
(C.21) 
(C.22) 
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+ sin Q!6 sin a-j cos as 
A 4 = 2(A 5 -Pa), 



\ 



1 - 



GT 



2t>Ai sin a e 



GT + HTv/JM 2 + v 2 sin 2 (3 



2v\ 2 sin 2 a-j 



A 6 = 2yAiA 5 cosa 6 , 
A 7 = 2A/A 2 A5Cosa7, 



(C.23) 
(C.24) 
(C.25) 



and the mass parameters according to Eqs. (|C.10|) -( |C.12|) . 

For later purposes, it is also useful to represent the parametrization in a slightly 
different form. Suppose now that A3, A4 are given parameters, in addition to Ai,A 2 . 
For A4, this is certainly consistent with the parametrization in Eq. ( |C.23| ) if A5 > 
and A 4 < 0. The former we assumed to be the case in order to get C violation, and 
the latter is always true in the MSSM. The constraint for A3, Eq. ( |U.22| ), then implies 
that A5 has a maximum allowed value for given A3, A4, 



< A 5 < A 



5, max 



2^3 + A4 



T 
G- 

v 



2vAiA 2 



cos a§ cos a.7 



— sin a 6 sin a 7 



\ 



1- 



GT 



GT + HTv/JM 2 + v 2 sin 2 /? 



(C.26) 



2vAi sin 2 a e J V 2t>A 2 sin 2 07 

This holds in the case that the expression on the right hand side is positive; otherwise 
no values of A5 are allowed (this typically happens for small v/T close to the minimum 
given by Eq. (|C.19|) ). The parameters A 6 , A 7 are still given by Eqs. ( |C.24j) , (|C.25|) . 



C.6. The effect of Higgs self-couplings and A t , ft 

In the analysis above, we ignored 1-loop effects from the SU(2) Higgses Hi, H 2 , and set 
At — $, — 0. Let us discuss here what happens when these assumptions are relaxed. 
Because we consider spontaneous C violation, A t , fi are assumed real. 

Clearly the introduction of A t ,fi,m 2 Si 7^ does not change the boundedness con- 



straints, Sec. p.2| . We will also not consider the condition of a local minimum, Sec. |C.4 



since this would be quite tedious. But looking for a stationary point as in Sec. |C.3| 
leads to useful observations. 

First, consider the effect of A t ,fi 7^ 0. Let us look at a local extremum constraint 
at some Vi,V2,0,<f>, obtained with mass parameters m\{T) , rn^iT) , m\ 2 {T) . By taking 
derivatives of Eq. ( |L2| ), we see that there is an extremum at the same point Vi, v 2 , 9, <j) 
also in the theory without any stop contribution in Eq. ( }4.2| ), but at the modified mass 
parameter values m\(T),m\(T),m\ 2 (T), where 



m\(T) = ml(T)- —h 2 t \jl\ 2 T [A + B cos 6 cos. 



,1/2 



(C.27) 
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ml(T) = m 2 2 {T) + ^-h 2 (l - \A t \ 2 )T [A + B cos 9 cos 0] 1/2 , (C.28) 

m\ 2 {T) = mUT) + — h 2 A t fiT [A + B cos 9 cos </>] 1/2 , (C.29) 

47T 



and A, -B are from Eq. 

We can now see that the values of, say, m\{T) + m\{T) leading to spontaneous C 
violation, differ typically from those obtained earlier on, rhl(T)+m 2 (T), by small effects 
~ -(3/(47r))/i t 2 (|i t | 2 + |/i| 2 )T[m 2 / (T) + (l/2)h 2 v^. (Recall that the dominant term 
in Eq. (|C.28|) which does not depend on A t , £i, was already included in our previous 
discussion.) Furthermore, the sign is negative, so that the part of the parameter space 
extending to the phenomenologically interesting region m\{T) + m\(T) > tends to 
decrease. The decrease can be rephrased by noting that A t ,ft ^ tend to decrease 
v/T, since they effectively decrease the coefficient of the cubic term which would be 
obtained from Eq. ( J4.2| ) in the limit mfj^) — > 0, and a smaller v/T makes C violation 
less likely. To summarise, we do not expect the inclusion of A t , pi ^ to change our 
conclusions. 

Similarly to Eqs. ( p.27|) - (|C.29| ), the 1-loop effects of the scalars Hi, H 2 are expected 
to change \m\(T) + m\{T^\jT 2 by terms parametrically of the type ~ X i mj/(2iTT). It 
is hard to dicuss this effect analytically, since in the general background of Eq. ( |3.4| ), 
the scalar mass matrix has the dimension 8x8. However, numerically we observe 
that the scalar contributions can also slightly increase the parameter space leading to 
spontaneous C violation, in contrast to A t , £i: parameters which would otherwise not 
result in a C broken minimum, can do so when the last term in Eq. ( |4.2j ) is included. 
Nevertheless, the effect is too small, numerically of order ~ Ajmj/(27rT) ~ 0.1, to have 
any qualitative significance. 



Appendix D. Integrating out the right-handed stop 



If the right-handed stop is heavy, it can be integrated out from the action in Eq. (2.1) 



In this case the electroweak phase transition is too weak for baryogenesis for physical 



Higgs masses in excess of 70... 80 GeV |14], [15], |16[ . Nevertheless, we summarise here 



how the couplings of the 3d SU(2) + two Higgs doublet model would change at 1-loop 
level if U is integrated out from Eq. (|2.1|), since we need the argument in Sec. 



5m\(T) = -A 7im[/ (T)T, (D.l) 

5m 2 2 (T) = -A 72m[/ (T)T, (D.2) 
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5m 2 



12CO 
S\ x 

5X 2 

S\ 3 
5X 4 
5X 5 
SX 6 
5X 7 



3 

4:71 
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712^,7 {T)T, 
T 



167r mu(T) 
3 T 



7?, 
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IQtt mu(T) 

3 T 

:7i72, 



87r mu(T) 

3 T 
87r mu{T) 

3 T 



1 7i2 1 



:7i 2 2 > 



167r mu(T) 

3 T 

:7i7i2, 



87r mu(T) 
3 T 
87r mu(T) 



72712- 



(D.3 
(D.4 

(D.5 

(D.6 

(D.7 

(D.8 

(D.9 

'D.10 



This integration is reliable in the limit that ^{1 ', XjjT ', g 2 T <^C mjj(T). As to the 
numerical magnitudes of the 1-loop corrections, recall from Sec. |2.1| that typically 



7i 



I712I "C 72 ~ 1. 
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